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Abstract 

Electrons with spin-orbit coupling moving in mesoscopic structures can often exhibit local spin 
polarization. In this paper, we study the influence of the Rashba coupling on the scattering of 
two-dimensional electrons from a circular disk. It is observed that spin-polarized regions exist, 
even if the incident electrons are unpolarized. In addition to the distributions of charge and spin 
current in the near-field region, we also analyze the symmetry and the differential cross-section of 
the scattering. 
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I. INTRODUCTION 



Spin-orbit interaction influences the electronic and transport properties of semiconduc- 
tors. For example, it lifts the degeneracy of the valence bands, modifies the electron g-factor,- 
and causes skew scattering in the presence of (spinless) impurities. Such a skew scattering 
is a possible mechanism for the extrinsic spin Hall effect.- In addition, spin-orbit interaction 
plays an important role in the recently proposed intrinsic spin Hall effect. 3 ' 4 It is also crucial 
in the mechanisms of spin relaxation and optical orientation in semiconductors.- 1 ^ 

It is highly desirable to generate flows of polarized spins in semiconductors with the 
help of the spin-orbit interaction. In these endeavors, the Rashba spin-orbit coupling^ in 
two-dimensional electron gas (2DEG) plays a special role since it allows manipulation of 
spin flows by varying the gate bias. This has motivated several creative proposals for its 
application.- To explore the possibilities, the effect of Rashba coupling in many types of 
mesoscopic structure have been investigated, such as a quantum wire,- a quantum ring,— and 
a quantum dot.— ^ In several studies, it was found that a device with a simple geometry, in 
combination with the Rashba coupling, could serve as a spin filter. For example, the device 
could be a T-shaped channel,— a quantum point contact,— parallel interfaces that cause 
double refraction,— or even just a curved wire.— After applying a magnetic field, we could 
further build a spin filter based on electron focusing,— or based on the an interferometer of 
the Stern-Gerlach type.— 

In this paper, we study the scattering of electrons by a disk in the 2DEG with Rashba 
coupling. 19 The analysis can be applied to a wide range of situations in which the radius R 
of the disk can be much smaller, roughly the same, or much larger than the electron wave 
length A. For example, the usual impurity scattering can be simulated with < A, while 
the scattering by an artificial mesoscopic disk corresponds to R ~ A. Here we focus on the 
latter case to search for possible effect of local spin generation. We find that, because of the 
Rashba coupling, spin-up and spin-down electrons are indeed separated and accumulate in 
regions of curved stripes. The associated charge and spin currents are analyzed in details. 
In addition, the properties of symmetry, as well as the differential cross-section, are also 
investigated.- 

This paper is organized as follows: Sec. II is the theoretical analysis. In Sec. Ill, the 
results from numerical calculations are presented. Sec. IV is the conclusion. 
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II. ANALYSIS OF THE DISK SCATTERING 



A. Hard disk scattering 

Consider a two-dimensional electron system with a circular disk at the origin, 

2 

where V(r) = Vq when r < R and otherwise. In the following, the potential of the disk is 
considered infinite ("hard" disk). The two-component wave function (ipi,ip2) T with energy 
E satisfies the coupled equations: 
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It follows that the energy eigenstates of the Hamiltonian with momentum k, helicity i] 
(E = % 2 k 2 /2m + rjak,rj = ±), and angular momentum (n+ 1/2) ft (where n is an integer) 
are 
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where can be a Bessel function, a Neumann function, or their linear combination, such 
as a Hankel function. We will choose the Hankel functions as the eigen-basis since their 
behavior at large radius suits the boundary condition for scattering. 

The energy eigenstate of the Schroedinger equation with energy E and angular momen- 
tum (n+ 1/2) ft can be written as 
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where E = h 2 k 2 /2m + ak = Ti 2 k' 2 /2m — ak'. The first two terms have positive helicity; while 
the other two terms have negative helicity. The most general eigenstate of the Schroedinger 
equation with energy E is a superposition of the ^ n 's, where the coefficients a n , b n , c n , and d n 
are determined by boundary conditions. Because of the circular symmetry of the potential, 



the angular momentum is conserved for each n-component. Therefore, each component can 
be considered independent during the scattering. 

The incident plane wave with momentum k and helicity rj can be decomposed as the 
following linear superposition: 
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The Bessel functions J n {kr) can be further decomposed as Hankel functions H^(kr) and 
H 2 (kr), the former correspond to outgoing circular waves, while the later correspond to 
incoming circular waves (with no phase shift). If an incident wave has a definite helicity, 
^in,+, then by comparing Eq. (4) with the components of Eq. (5) at large distance, one will 
obtain b n = i n /2y/2 and d n = (i.e. no incoming circular wave with negative helicity) for all 
n. The coefficients a n and c n need to be determined from the boundary condition at r = R. 
For a hard disk, it can be shown that 
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where k = kR and k = k'R. Notice that the nonzero probability amplitudes c n lead to 
outgoing waves with flipped helicity. For reference, if the incident wave is ^m,-, then b n = 
and d n = i n /2\^2. At the mean time, the roles of k and k', as well as the roles of a n and c n , 
have to be interchanged. 

A note on the unitary condition: for convenience of discussion, consider an incoming wave 
with positive helicity and angular momentum (n + l/2)h (i.e., the 6 n -wavelet). Because the 
angular momentum is conserved during the scattering, the electron can only be scattered to 
a n and c n channels with the same n. From particle conservation at large distance, one expects 
that the probability amplitudes in Eq. (4) should satisfy \a n \ 2 + (k/k')\c n \ 2 = \b n \ 2 (= 1/8) 
for all n, which has indeed been confirmed in our numerical calculation. 

B. Properties of symmetry 

The system has a mirror symmetry with respect to the x-axis. Therefore, by analyzing the 
Schrodinger equation with y replaced by — y , one finds ^(r*) = —a y ^/(f), where r— (x,y), 



and r* = (x, —y) is the mirror-reflected point of f. Such as relation can also be obtained by 
a space inversion of the (three-dimensional) coordinate, followed by a rotation with respect 
to the new y-axis by 180 degrees. Consequently, for the expectation value of the spin, we 
have 

(5,(0, Sy(n, S z (n) = (-S.(r), S y (r), -S z (r)). (7) 
In a Rashba system, the current density operator is defined as, 

j = -f --^xR (8) 
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Therefore, the distribution of the expectation value of the current density has the following 
symmetry: 

Ux(r),j y (f)) = (j x (f), -j y (r)). (9) 

We adopt the generally accepted definition of the spin current density operator,— j] = 
whose expectation values have the symmetries, 

u*(njy(n) = (-mjm- ( io ) 

These symmetries will be confirmed by the numerical results in Sec. III. 



C. Asymptotic behavior of the scattered wave 



For convenience, the wave function can be seperated into an incident plane wave 

and a scattered wave. At large distance with kr 1, the scattered wave has the asymptotic 
form, 
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It can be shown that ft • g = 0. Also, for the incoming plane wave ^m,+, 
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^2fi#fi = |f| 2 (sin0, -cos0,O), 
i=i 
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J29i°9i = |g| 2 (-sin0,cos0,O). (13) 
i=i 

Therefore, f-spinor and g-spinor possess spins with opposite directions at large distance. 
Both spins lie on the plane and are perpendicular to the direction of propagation. One can 
obtain the same equations for the incoming wave ^m,-, but the signs of the spin expectation 
values are opposite. 

After a straightforward calculation, one can show that the scattered current density at 
large distance is 

L = -(- + j)\t\f + -(--j)\t\f, (14) 

r \m h J r \ m n J 

from which the differential cross-section a' ((f)) = r\j sc \/\jm\ can be calculated. For incoming 
waves \l/m,+ and ^m,-, the current densities \jm,+ \ and \jm,-\ are hk/m+a/h and Tik' /m—a/Ti 
respectively. In fact, they are equal in magnitude if the two incident waves have the same 
energy. Therefore, the differential cross-sections a n for incoming waves with helicity rj are 
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where f v is the f-spinor in Eq. (12), but with the coefficients a n and b n suitably chosen for 
the scattering of ^m,^, similarly for g,, (see the discussion following Eq. (6)). The differen- 
tial cross-sections a' and cr 1 represent heli city-preserved and helicity-flipped scatterings 
respectively. If the incoming wave is an incoherent mixture of both helicities with fractional 
populations P v , then the differential cross-section is simply the weighted average of the two 
differential cross-sections: a' = P+&+ + P~o'_. 



III. DISTRIBUTIONS OF CHARGE AND SPIN 



In the following, we report on the distributions of spin density, charge current density, 
and spin current density, assuming the wave length of the incident wave A = R. We have 
also studied the cases with a larger A (e.g. A = 3R) and a smaller A (e.g. A = R/3). These 
results are not presented since the main difference is the change of scales. In the limit of 
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FIG. 1: Distribution of out-of-plane spin S z , with a = (l/10)h 2 k/2m. The incident wave is an 
equal mixture of both helicities rj = ±. Lighter (darker) regions represent the regions of spin up 
(down). The peak intensity in this figure is approximately 0.22 (1 for fully polarized). 

A ^> R, whose scale is more relevant to the case of impurity scattering, only the components 
with the smallest angular momentum (n=0 and —1) need to be considered. From Eq. (jSj), 
one finds that both |co| and |c_i| —>■ as kR <C 1. Therefore, there would be little change 
of helicity in the long wave length limit. 

For comparison with realistic values, we choose m = 0.068m e for electrons in the GaAs- 
AlGaAs heterojunctions. The radius of the disk is fixed at i?=1000 A. The corresponding 
Fermi energy and electron density for Xp = 1000A are 2.2 meV and 6.3 x 10 10 /cm 2 respec- 
tively, which are typical values. To enhance the visual effect of the spin-orbit coupling, the 
Rashba energy ak is chosen to be one-tenth of the kinetic energy h 2 k 2 /2m, which requires 
a = 0.35 nm-eV, about one order of magnitude larger than the value in GaAs. 

The out-of-plane (z) spin component results from a spin-unpolarized incident wave with 
an equal (incoherent) mixture of both helicities is plotted in Fig. 1, which is antisymmetric 
with respect to the horizontal x-axis: S z (f*) = —S z (r). Notice that the incident waves with 
opposite helicities (but the same energy) have different wave vectors (A/ — k — 2ma/h 2 ). 
Therefore, their interference patterns for opposite helicities with spins point at opposite 
directions are slightly displaced with respect to each other. Because of such a displacement 
between r\ = + and r\ = — , regions with net z-spin still exist after partial cancellation. The 
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FIG. 2: Distribution of charge current density near the disk. The incident wave is ^i n ,+ with 
wavelength A = R. [a = (l/W)h 2 k/2m\ 

existence of S z relies on the scattered part of the wave function in the near-field region [see 
Eqs. ( ITT1) and ({TBI ] and would decay to zero at large distance. 

Unlike the spins in Fig. 1, the distribution of charge current (Fig. 2), as well as the 
differential cross-section (Fig. 3), which is defined using the charge current densities, are not 
sensitive to the strength of the Rashba coupling, and look very similar even if the coupling 
is turned off. Notice that Rashba spin-orbit coupling in fact preserves the helicity. The 
major cause of the helicity flip is the potential V(r), which is incompatible with the helicity 
operator. In Fig. 2, the distribution of charge current density from the scattering of ^i n ,+ 
shows the expected pattern of the flow around the disk. If the incident wave is a mixed state, 
then instead of cancellation, the slightly displaced current densities from both helicities will 
add up. 

The differential cross-sections for \l/i n) + are shown in Fig. 3. It can be seen that the 
helicity-preserved scattering {o~' ++ ) peaks at the forward direction (0 = 0,2tt), while the 
helicity-flipped scattering (er+_) peaks at the backward direction (jp — tt). At the backward 
direction, a' ++ {it) = 0. Therefore, the helicity of the electron has to be flipped, but its spin 
remains conserved. At longer wave length A = 10R, cx' + _ gains more weight and the total 
differential cross-section o~' + = a', , + becomes more isotropic. 

The distributions of spin current density are shown in Fig. 4.— A prominent feature in the 
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FIG. 3: (Color online) Differential cross-sections (in units of R) for wave length A = R (solid line 

for a' + , , dashed line for cr' | ) and A = lOi? (dash-dotted line for &' ++ , dotted line for cr' ). The 

numbers in the legend are total cross-sections after integration over ang le. [a = (l/W)h 2 k/2m\ 

figures is the overall trend for the j%(r) vectors to point to the left. This counter- intuitive 
behavior is simply related to the fact that the spin current density is equal to the product 
of velocity and spin, where the velocity is in the positive x-direction and the spin points to 
the minus ^-direction for \I/in,+, therefore the scattering j y s vectors generically point to the 
left. On the contrary, for with spin points to the positive y-direction, the direction 

of the flow will be reversed. This also explains why the magnitudes of the spin currents 
jx,z are sma ]i j n m ost of the regions, since the original incident current has no S x and S z 
components. 

In Fig. 4, it can be seen that j* oscillates both in amplitude and direction between the 
curved stripes. An incident wave with opposite helicity would reverse such a flow. Therefore, 
part of these flows are cancelled if the incident current is not polarized. However, local spin 
current that oscillates in space still exist, similar to the case of the spin density in Fig. 1. 
All of the local spin currents ' y,z would vanish for unpolarized incident electrons once the 
Rashba coupling is turned off. However, if a is nonzero but the disk is removed (i.e., free 
space), then there exist spin currents j% = a/2x, j v s = —a/2y, and j z s = 0, which is the 
background spin current cautioned by Rashba.— It reflects the unsatisfying current status 
on a proper definition of the spin current. 
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Our system with unilateral current flow seems to be the same as the 2DEG driven by an 
electric field,- but there is no global transverse spin current j*, no matter the disk is removed 
or not. This does not contradict the result of the proposed intrinsic spin Hall effect in a clean 
2DEG. 4 The incident charge current flowing to the right can be understood as originating 
from the slightly unbalanced electrochemical potentials on the two leads far away. In our 
case, all the electrons are moving along the direction of the potential gradient, instead of 
moving at all directions on the Fermi surface in Sinova et aPs paper.- Consequently, no spin 
Hall effect is expected if one follows similar semiclassical analysis in Ref. |4. 



IV. CONCLUSION 



In summary, the influence of the spin-orbit coupling on the scattering of 2D electrons 
from a hard disk is studied. Such a simple setup offers us a good opportunity to investigate 
the properties of spin and spin current in details. We focus our attention on the near-field 
regime, where the scattered wave is comparable to the incident wave, and appreciable out- 
of-plane spins can be found. This work offers us a clear understanding of the microscopic 
dynamics around the mesoscopic disk, and could serve as a basis for future works considering 
a threaded magnetic flux in the disk, or a hybrid device involving a disk as a component. 
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